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Abstract 

Starting from the dipole representation of small- a; evolution we implement the running of 
the coupling in a self-consistent way. This results in an evolution equation for the dipole density 
in Borel (b) space. We show that the Borel image of the dipole density is analytic in the 
neighbourhood of b = and that it is equal to the BFKL solution at b = 0. We study the Borel 
singularity structure of the dipole cascade emanating from a virtual photon at small x and find 
a branch cut on the positive &-semiaxis starting at b = l//3o- This indicates the presence of 
1/Q 2 power corrections to the small- a; structure functions. Finally we present numerical results 
in the context of D.I. S. 



PACS numbers: 11. 55. J, 12.40.M 



1 Introduction 



The study of QCD in the high energy limit has a history spanning at least twenty years and is still 
actively pursued. Within perturbation theory (pQCD) an understanding has emerged concerning 
hadronic processes in the semihard regime, such as near-forward elastic and diffractive scattering, 
small- x deeply inelastic scattering and rapidity gap events. The evolution of the partonic cas- 
cade stemming from the initial state hadrons is known to be described in the leading logarithmic 
approximation LLA(x) by the BFKL equation either in its original form or in its improved 
reformulations However, the picture is yet incomplete because of two major problems encoun- 
tered in this regime. The first has to do with unitarity ||, which is violated by the LLA(x) result. 
The second concerns the infrared sensitivity of the LLA(x) result, which is infrared finite but does 
receive contributions from low momentum regions, where physical observables are sensitive to large 
non-perturbative corrections. It is the second problem that we shall consider in this paper. 

Unlike the more conventional high momentum transfer {Q) processes with final state particles in 
the central rapidity region, the semihard kinematic regime is characterised by three ordered scales, 
y/s >Q> Aqcd- Note that the use of pQCD is justified by the second of the previous inequalities, 
in the absence of which, Regge phenomenology of soft processes is the only available approach so far. 
The pQCD factorisation theorem, applicable for high-Q processes, states that inclusive observables 
can be expressed as a product of parton distribution functions and a hard scattering cross section 
(matrix elements and coefficient function in OPE parlance). In the semihard regime this theorem 
needs amendment. Here, the product is not only a convolution in longitudinal momentum fractions 
of the partons but also in their transverse momentum. This is known as fc^-factorisation and 
in impact parameter space it translates to a convolution in impact parameter of the radiated 
secondary partons and an off-shell hard amplitude or cross section. This leads us to the first of 
the two theoretical inputs of the present study, namely the use of the impact parameter or dipole 
emission formulation of semihard processes, as developed by Mueller H and independently by 
Nikolaev and Zakharov ||. The merit of this formalism, apart from the simplicity of the final 
results, is that it organises the process in terms of sequential soft gluon emissions in the s-channel, 
which in the Coulomb gauge admit a clear physical interpretation. They are to be understood as 
initial state radiation occurring before the hard scattering. Final state interactions are not included 
in the evolution equation because they only lead to a unitarity rearrangement that does not affect 
inclusive observables. 

To study the infrared sensitivity of a semihard process we are bound to consider a subset of next- 
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to-leading logarithmic corrections, specifically the ones that have to do with the running coupling. 
In the BFKL approach, the evolution kernel, although infrared finite, is sensitive to the running 
coupling because of diffusion towards the infrared of the transverse momentum along the exchanged 
gluon ladder. With specific assumptions about the infrared dynamics (i.e. long wavelength gluon 
propagation) and the scale of the running coupling, the effect on the Regge trajectories has been 
analyzed in refs. Q. In the present study we take the point of view that in the dipole formalism, 
just like in timelike parton cascades, the scale of the running coupling is set by the virtuality of the 
emitted gluon. The problem of translating this into impact parameter space is resolved through 
the use of the Borel transformation. We require that the dipole evolution kernel has the same 
Borel singularity structure as the one obtained in momentum space. This condition is sufficient to 
determine the scale of the coupling as a function of the impact parameters of the emitted gluon. 
Therefore, the second theoretical input in our approach is renormalon analysis Q. Once the dipole 
kernel with running coupling is constructed, we study the power corrections that are generated 
when this kernel is used to evolve the wave function of an initial hadronic state such as a virtual 
photon in small- x D.I.S. 

The structure of this paper is as follows. In section 2 we highlight the main concepts and results 
of the dipole formalism. Here we also review the derivation of the dipole evolution equation with 
running coupling. In section 3 we discuss the Borel singularity structure of the dipole density. We 
find a series of infrared renormalon singularities at b/3o = n, n = 1,2,3,... as well as a leading 
singularity at bfio = 7, with < 7 < 1 the anomalous dimension of the solution. We show that the 
Borel image of the dipole kernel is analytic in the neighbourhood of b = and compute its action 
as a power series expansion in b. In the limit 6^0 the BFKL result is recovered. Iterations of 
the kernel turn the renormalon poles into branch cuts but do not shift their positions. In section 
4 we discuss small- x D.I.S. in the dipole formalism using the results derived previously. The main 
conclusion is that the leading Borel singularity is at b = l/f3o, which implies 1/Q 2 power corrections 
consistent with the Wilson OPE expectation. The dipole evolution equation for fixed coupling can 
be solved by virtue of its scale invariance and leads to the same spectrum as that of the BFKL 
equation. For the Borel image this scale invariance is necessarily lost and the resulting evolution 
equation can only be solved in numerical approximation by successive iterations. In section 5 
we present results for the first few iterations of the Borel transformed kernel in deeply inelastic 
scattering, and demonstrate that the singularity conforms with that discussed in section 4. We 
summarise in the final section. 
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2 The dipole kernel with running coupling 



We consider small- x D.I.S. as a specific example of a semihard process and describe it in the rest 
frame of the nucleon target of mass mjv- Let (z, r) denote the wave function for 7* to fluctuate 
into a q-q dipole of transverse size r and with the quark carrying longitudinal momentum fraction 
z. Then, the virtual photoabsorption cross section a(x, Q 2 ) for given polarisation of 7* (T, L), can 
be written as []|, ||] 

<?t,l{x,Q 2 ) = J dz J d 2 r\^ L (z,r)\ 2 a dN (Y,r), (1) 

where a^N denotes the dipole-nucleon total cross section and Y = ln(z/x) is the large rapidity 
parameter to be resummed within pQCD. The physical interpretation of eq. (Q) is that the tran- 
sition 7* — * qq, with light cone time scale r = 0(l/xmisr), occurs much earlier than subsequent 
interactions with the nucleon target. The lowest order in a s dipole cross section is 

*2k*) = ^ / ^ 2 V(l)Re(l ~ e" iIr ) , (2) 

with V(l) proportional to the absorptive part of the gluon-gluon-nucleon-nucleon vertex function. 

The LLA(x) radiative corrections are generated by emission of soft gluons by the initial q-q 
dipole with longitudinal momentum fractions strongly ordered as 

z > z\ > z 2 » ... > z n . (3) 

The light cone time scale for the fc-th gluon emission is 17 = 0(^Q/k?), therefore it is emitted 
independently of the previous ones (r^ 3> r^+i). The 7* wave function ip( n \z, Zi,r,ri) for the 
emission of n soft gluons with longitudinal momentum fractions Zi and at transverse positions 
can be calculated order by order in pQCD and the inclusive probability distribution for 7* to 
fluctuate into q-q plus soft radiation, <£(z,r), is obtained by the square of integrated over 
the phase space of the emitted gluons and summed and averaged over their polarisations {Aj} and 
colours {aj}. 

n j 2 r*'—x A 

* (n W) = E U^r r ^i^r 1 ^^,-^^,^,...^)! 2 , (4) 

\i,a,i i=l 1 

&°\z,r) = |^ (0) (^,r)| 2 . 

Here and until section 4 we suppress the 7* polarisation indices (T, L) for brevity. Indeed it is the 
projection of these polarisation states for the virtual photon that is responsible for the fact that 
$(°)(z,r) is taken to be a function of r only. 
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In the Coulomb gauge the soft radiation can be viewed as a cascade of colour dipoles emanating 
from the initial q-q dipole. It is useful to introduce the dipole density n(Y, r, p) [||, [lC| 1, for emission 
from the initial q-q dipole of size r of a dipole of size p and with the smallest longitudinal momentum 
fraction in the emitted dipole bounded from below by e~ Y . The dipole density is independent of 
the projectile bound state properties. When convoluted with the 7* wave function squared it 
generates the number density N(Y, p) of dipoles of size p and rapidity larger than Y as 

N(Y, p) = J dz J d 2 r $ (0) (z, r) n(Y, r, p) . (5) 

The whole dipole cascade can be constructed from the repeated action of a kernel JC on the initial 
density no(Y,r,p) through the dipole evolution equation 

pY poo 

n(Y,r, p) = n (r, p) + / dy dr'/C(r, r') n(y, r, p) . (6) 



with driving term 

no(r, p) = n(Y = 0, r, p) = r5(r — p) . (7) 

The evolution kernel /C can be calculated from the one soft gluon emission probability & l \z, r). 
The result is 

•« ( ,, P) = *<°> M fff^ag.^, (8 , 

where hatted vectors will always denote the distance from the second parent emitter, f 1 := ri — r. 
The full kernel /C is 

JC{r,r) = —S{r-r) dr —„ —3 + 2—^ / dr J(r, r , r )-=- 5 . (9) 





The first term of K, comes from the one loop virtual corrections to the initial q-q dipole of size 
r. The 5-function denotes the absence of gluon radiated into the final state. The second term 
corresponds to real emission and can be read off from eq. (||) taking into account eq. (|6|). A factor 
of 2 is included here to account for the fact that after one dipole splits into two, each one of the 
offsprings can act as parent for subsequent emissions. The change of integration variables 

dV = J(r,r',f')dr' dr (10) 

generates the (triangle) Jacobian § 

poo 

J(r, r', f') = 27T r' f' / dun Jo{Kr) Jo(rer') Jq{kt') . (11) 
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Unlike the BFKL kernel, both virtual and real parts of the dipole kernel K, are separately IR 
finite in the limit r',r' — > oo. For fixed a s and in the limit N c — ► oo, for which 2Cf — ► N c , K, is 
conformally invariant and has the same spectrum as the BFKL kernel, i.e. 



dr%(r,r')(r' 2 y = ^£ x ( 7 )(r 2 )^ 



7T 



where 



x(7) = 2*(l)-*(7)-*(l-7) 



(12) 



(13) 



is the BFKL spectral function. Indeed, the two approaches lead to the same phenomenological 
results for inclusive observables. The equivalence between the dipole and BFKL formalisms at the 
level of light cone perturbation theory diagrams has been shown explicitly in ref. plfl. 



After the introduction of the dipole density n, the dipole cross section GdN of eq. (^) takes the 
form 

d 2 p 



a dN (Y,r) 



n(y,r, p) crJJ(p) . 



(14) 



Eqs. ([l], |l4| ) provide the factorised form of small- x D.I.S. in impact parameter space. We note 
that dependence on target mass scale enters through <r^ and that the universal part of the 
photoabsorption cross section a(x, Q 2 ) is indeed the dipole density n, which is independent of 
projectile or target. 

The problem of how to include the running coupling in the dipole evolution kernel /C(r, r') has 



been studied in ref. [12]. The form of the kernel with running coupling is fixed by the following 
two requirements. The first is that in the emission probability of a soft gluon with transverse 
momentum k the coupling runs as a s (k 2 ). The second is that the emission probability in impact 
parameter space should generate the same singularities in Borel space as the Borel transformed 
emission probability in momentum space. Then the form of the one-dipole emission probability 
<&( 1 )(£,r) is completely determined. Explicitly it is 

* {1 W) = ^ (0 W) f— /d 2 n rdrdfJi(T)Ji(f) f^j^- > 

vr d J z\ j Jo Jo ur /i: (l — uj) 



1/2 



where 




A 2 (r,f, W ) = (r 2 r(f 2 )( 1 -), 



r 2 r{f 2 )^ 



(15) 



(16) 
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In the limit of fixed a s the above equation reproduces the expression for the one-dipole emission 
ffrW of eq. (||). It turns out that in impact parameter space the scale of the coupling is not a simple 
function of the distances n, r\ but is weighted with three further parameters r, f and u. 

The Borel image ^^'(z, r;b) is defined, as usual, by the transformation 



$W(z,r) = / db& 1 Hz,r;b)e- b / a °( Q2 K 
Jo 



(17) 



This transformation can be easily inverted when the one-loop running coupling is used. Then the 
Borel image a s (k 2 /Q 2 ;b) of a s (k 2 ) is 

-bfo 



k_ 2 

Q 2 



a s (-j;6) 



Q 2 



(3 = — (—N c --Nti. 



(18) 



and from eqs. (15,17) we obtain 

* (1) (2,r;6) = ^<I> (0) (z,r)ln(z) / d 2 n /°° drdf Ji(t) Ji(f ) /' 1/2 U/ 
vr d 7 jo Jo tj 1 /^ (1 - a;) 1 / 



r 2 f 2 



2 -2 



rr r 



l ' l 



2 + (f 2 ) 6/3o r 2 + 2ri . ?1 ( r 2 ) 



2w6ft)/ f 2\(l-w)fe/3 



(19) 



The integrals over 2i,ri, r, f and a; can be performed and the result is 



vr w r(l + 6/? ) 



(20) 



From eqs. ([j^, 20) we can construct the dipole evolution kernel in Borel space JC(r, r'; b). It consists 
of two pieces. The virtual contribution is the coefficient of the 3>(°) in eq. (|20|) with the sign inverted 
as required by unitarity j|. The real contribution is obtained by the coefficient of in eq. (|i~9|). 
After renaming r\ — ► r', changing integration variables as in eq. ( JTo| ) and performing the integration 
over r' we obtain the explicit form of the kernel. It is 1 



N c (Q 2 r 2 \ m) 



K(r,r';b) = ~r . , 

7T \ 4 



+ 2 



iV c 1 Q 2 r 



I r(l + 6/3 ) 

2„,/2 \ rl 



7T 2 r' 



„/2 



5(r - /) 

da; r(l -Lobf3o)T(l - (1 -u)b/3p) 
o ^^(l - u) 1 / 2 T(l + w&A,) T(l + (1 - u )bPo) 



b(3 -l 



2 F 1 1 - 6/Sb, 1 - 6A); 1 



' -r.2 



x The factor of 2 in front of the second term of this equation has been inadvertently omitted in eq. (16) of ref. |lS 
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/ r 2_ r /2\ / 2 \ / 2 \ 

( 2 \ i^6/3o / 2X1 

Plj 2^i^6/3o,-^o;i;^J| , (21) 

where r< = min(r, r') and r> = max(r, r'). Finally, the evolution equation for the dipole density 
in Borel space reads 

—n(Y, r, p;b)= I dr' \ dbX(r, r'; b') h{Y, r' , p; b - b') , (22) 
with boundary condition 



o Jo 



n(Y = 0,r,p;b)=r6(r-p)S(b). (23) 



3 The Borel singularity structure of the dipole cascade 



To study the Borel singularity structure generated by the dipole kernel K, we consider its action on 
a test function of the form (r 2 ) 7 . Defining the function x(7> b) by 

b/3 



7T 



(24) 



we obtain the expression of xili b) from eq. (|2l|), 

T(-bpo) 



X(l,b) 



+ 



r(i + bpo) 
r(- 7 -6A)) i 

T(l + 7 + bPo) vr 7 w i/2(i _ W )V2 r(l + w6/3 ) r(i + (1 - 



da; 



r(i - a;6/3 ) r(i - (l - u)bp ) 



r(i + 7 ) r(6/? ) 



r(l + 7+ (1 -£«;)&#,) T(l + 



(25) 



L r(- 7 ) r(i - 6/3 ) r(i - 7 - (i - w )&/3b) r(i - 

where the first term comes from virtual correction and the second from real emission. The virtual 
contribution to x(l>b) contains a series of poles at b/3o = 1, 2, 3... which are identified with the 
IR renormalons and correspond to power corrections of 0{(m? N /Q 2 ) n ), n = 1, 2, ... Note that these 
poles are independent of the specific form of the test function. This set of poles, resulting from 
the exponentiation of soft radiation, has also been derived in the context of the Drell-Yan process 
in ref. p3|. In addition, we observe the presence of a series of poles at 6/?o = n — 7, n = 0, 1, 2... 
generated by the r(— 7 — 6/?o) dependence of the real contribution to xil-, b)- For i?e( 7 ) > m these 
poles correspond to IR renormalons for n > m. Their IR origin is established by observing that 
these singularities arise from the r' > r integration region of eq. (|2"2|), where the offspring dipole is 
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emitted with size larger than the parent dipole. Taken at face value the 7-dependent poles indicate 
the presence of 0{{m 2 N / Q 2 ) n ~^) power corrections. For the moment it suffices to note that in this 
section we are studying the dipole cascade independently of its embedding in a particular physical 
process. The interpretation of the 7-dependent power corrections in terms of OPE matrix elements 
is not an issue here. 

Inspecting the evolution equation (^) we now examine whether the kernel introduces singularity 
at the lower end of the ^'-integration. The answer is that for fixed r the limit b — > of /C(r, r', b) is 
not singular. The 6 = pole appearing in the virtual contribution to JC is of UV origin and cancels 
against a corresponding singularity in the real contribution. The presence of such a singularity 
can be seen already in the fixed coupling case. The integral in eq. (^) is logarithmically divergent 
in the region 77 — > 0. One way of tracing the cancellation of the UV singularities is to impose a 
lower cutoff in the impact parameter integration as in ref. Q]. Alternatively, in the case of running 

coupling we can use the Borel parameter 6 as a dimensional regulator. Defining 

r 2 

C = ~f (26) 

r > 

we note that the UV divergences in the real part of K, arise from the region £ — ► 1. This is the region 
where the gluon is emitted at impact parameter almost equal to that of the parent dipole. In the 
limit 6^0 the only UV divergent contribution comes from the first term inside the curly brackets 
of eq. (^l]). It is due to the divergence of the hypergeometric function 2-^1(1 — &Ab 1 — bfio; 1;£ = 1) 
at b — > 0. To parametrise this divergence we use the linear transformation [14] 

a Fi(l - bp , 1 - bp ; 1; = (1 - C)- 1+2b ^ 2 F 1 (b/3 , bfo; 1; • (27) 

Then, in the limit b — > the kernel takes the form 

jC(r,r';b) = ^ ( ° \-±—5{r - r') 




(l-C)- i ™ 2 F 1 (0,0;l;O + (UVregular)j . (28) 

The singular factor (1 — Q~ 1 + 2b Po can fog expressed in terms of distributions as 

1 



(1 _ c) _i +aW) = 1 §{1 _ c) + 
zopo 



i-C 



+ ... (29) 



Substituting this into eq. (28) it follows that the UV singular term in the virtual part of tC cancels 
against the one in the real part. Since the kernel is regular around b = we compute the power 
series expansion in b of x(7, 6). From eq. (25) we obtain 

X(7, b) = X (7) + b(3 oX {1) (7) + 0{b 2 dl) . (30) 
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where 



X (1) (7) = --X(7) - 2*(l)x(7) + \x{l? + \x{l) (31) 
7 11 

The 0(l/b) singular terms have cancelled as anticipated and the O(b ) term is the BFKL spectral 



function x(l)i ec l- (13). 



Returning to the dipole evolution equation (22) we introduce the usual Laplace transform with 
respect to the rapidity Y 

nUr,p;b)= dYe-" Y n(Y,r,p;b), (32) 



and the anomalous dimension 7 via the Mellin transform 

f°° dr 2 ( r 2 \ 7 

n un (b) = —2-[~2 I ^u>(r,p;b). (33) 



Then eq. (22) with the boundary condition ( [23] ) through the use of (24) becomes an equation for 
the spectral amplitudes h ua (b), 



l2 „2 \ h 'P° 



un^(b) = 25(b) + / db'xil - b%, b') p/- ] ^ >7 _ 6 ^ (6 - b') . (34) 



7T Jo 

This equation is universal for all semihard processes. The scale Q 2 is only constrained by the 
requirement Q 2 3> Aq CD and its specific value is otherwise arbitrary. It is identified with the 
7* virtuality here because we consider the specific case of small-x D.I.S. The dipole density is 
reconstructed from the spectral amplitudes as 



(35) 



Eq. ( pip is an integral equation of the Volterra type. Such equations are known not to have 
eigenfunctions for bounded kernels. One way of constructing solutions is by iteration. In our case, 
the iterative solution is given by the formal expression 



00 

»M = £ ng* (6) , (36) 

k=0 



with 



and 



n£» = nW(b) = I {^f \ ° xfr-bM, (37) 

2 / N \ k /0 2 n 2 \ b ^° r b 
^ {b) = Jo d6 fc-ix(7-(6-6fc-i)A,b-6fc-i) 
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x(7 ~ (b- &fc-2)A), h-i ~ h-2) ■■■ 
db 2 xil ~ (b- b 2 )(3 , 63 - b 2 ) 
62 d6i x(7 - (6 - &2 - 61) x(7 - ftA), &i) , (38) 



for k > 2. The index counts the number of iterations of the kernel in eq. (34). 



Having identified the IR renormalon poles of x(7> b), eq. (|25|), we can study the Borel singular- 
ities of the spectral amplitude n U) ^(b) from eq. (^). The contribution to n^)y(b) from the virtual 
(7-independent) terms only is 

2 (N c \ k (Q*f?\ bfh /■*> r(-(6-6 fc _ a )A)) Z^- r(-(6 fc _! - 6 fc _ 2 )/3 ) 



V 4 / ^ "- i r(i + (6-6 fc _ 1 )A))yo " r(i + (& Jfe _i + & fc _ 2 )A ) r" 
fb2 J , T(-(6 2 - 60/3Q) r(-6^o) 

5o 



1 r(i + (6 2 + 6 1 )/? )r(i + w ' 



Because the 6-integrals are nested, for b < l//3o no new singularity is introduced in the solution from 
these virtual corrections. For b > 1//3q some of the T- function poles will be encountered on the b- 
integration paths turning the leading b = 1/(3q pole into a branch cut. Real contributions introduce 
singularities on the positive 6-semiaxis for real 7, and for < 7 < 1, the leading one is found from 
eqs. (37, [3|) to be at b = 7//?o- We note that for bj3o < 7 the real contributions to n^(6) will not 
generate singularities because none of the poles of the x-functions in eq. (|38|) are encountered along 
the nested 6-integration paths. For bfio > 7, as in the case of virtual contributions, poles are found 
along some of the 6-paths turning the leading b = j/(3o pole into a branch cut. From eqs. (^, |38| ) 
it follows that there are also singularities for negative values of b. This means that b = 7/^0 is 
the position of the leading singularity for positive b only when 7 is in the range < 7 < 1. If we 
were to take 7 in the range n < 7 < n + 1, n £ Z, then the leading IR singularity would be at 
7 — n, and the singularities for negative b would be displaced further to the left on the negative 
6-semiaxis. Although these singularities do not affect the estimate of the power corrections from 
the non-perturbative effects, they reflect the fact that as 7 increases the dipole densities generated 
from these test functions become increasingly IR divergent. 

The above argument can be repeated for the contributions to the solution from cross products of 
virtual terms for some of the x-functions and real terms for the rest. The net result of this analysis 
is that the region of analyticity in Borel space of n wri {b) with < 7 < 1 contains the interval 
< b < 7//?o and the leading singularity is a branch cut at b = 7//?o- I n the asymptotic limit 
Y — > 00 the anomalous dimension is 7 = 1/2 + iv with v the spectral parameter to be integrated 
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over in eq. (p5|). Such a Re(j) implies the presence of power corrections of 0{mfq/Q) for the 
dipole density n(Y,r,p). These 1/Q power corrections have also been reported in the context of 
conventional BFKL approach in ref. Jig]. 



The evolution equation (p5) is not conformally invariant because the action of the dipole kernel 



results in a shift of 7 by b' (3q in eq. (34). However, having established above a region of analyticity 
in Borel space that contains b = we can use the 6-expansion as a measure of deviation from the 
conformal limit. From eqs. |38| ) the 6-expansion of the spectral amplitude generated after k 
iterations of the kernel is 



2 fN r 



b k- 



uj \ttuj J I {k — 1)! 

b k (k 



x(i) k + 



(k-iy. 



A) In 



X(7) fe + /?0X (1) (7)X(7) 



k-l 



2 (k- 1)! 



U /3ox(7) fc - 1 ^ + 0(6 fe+1 ) 



d7 



(40) 



The full spectral amplitude n^^{b), defined in eq. (36), can be constructed in the small-6 region 
from summing the powers of b. Summation of the leading powers of b, i.e. retaining only the 
0(6 fc_1 ) terms in eq. ( |40|) is required in the region 



6< 1, 



1. 



TTUJ 



(41) 



This 'resummation' of the 0(N c x(j)b/ '(ujtt)) terms is well defined if the 0(b k ) terms are small 
relative to the C(6 fc_1 ) terms in eq. (|40|). For the 0(b k ln(Q 2 p 2 /4)) term this translates to the 
requirement 

CVr? / 1 1 \ 

< 1 • (42) 



6A) In ^r- < 1 6 



,a s (Q 2 ) a s (l/p 2 ), 

Hence, leading 6-power resummation is valid in the conformal limit of fixed coupling. In this case 
the Borel transform can be readily inverted in terms of a fixed a s to yield the well known answer 
for the BFKL spectral amplitude 



(JXJ r) 

E- 



1 N c a sX {l) 



ITU 



lTLd 



+ ... 



lo - N c a s x{l)/n ' 

and via eq. (p5|) the asymptotic solution for the dipole density is obtained (l( 



(43) 



n(r, p,Y) 



with ap — 1 = A.(a. s N c /n) In 2. 



1 r exp[(ap — l)y] 
2p V(7/2)«,iV c C(3)y 



exp 



7rln 2 (r//j) 
14a s iV c C(3)y 



(44) 
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4 Small- x D.I.S. and power corrections 



The structure functions in the small- x region can be written in the following form 

F T ,L^,Q 2 ) = IT— t dz I d 2 r^ L (z,r)a d:N (Y = ln(z/x),r). (45) 
Using the definition of the dipole-nucleon cross section, eq. fll4|), the above expression becomes 

FtA^Q 2 ) = ^-£dzJd 2 T^(z,r)J^n(X,r,p)a ip,m N ) 

Q ~ 1 d2p ^N^ L (Y,p)a (p,m N ). (46) 



47ra em J 2np 2 

This equation is the factorisation theorem for the structure function in the small- x regime. As 
mentioned in the introduction, it contains in addition to convolution in the longitudinal momentum 
fraction a convolution in impact parameter. <&(°) is the lowest order transition probability for 



7* — ► qq. (This is the lowest order impact factor of ref. [17|.) It can be calculated reliably in 
perturbation theory in the limit Q 2 3> Aq CD . The dipole density n(Y,r,p) contains the universal 
evolution of the initial q-q dipole into a cascade of dipoles with ordered rapidities. This is also 
calculable in perturbation theory as we saw in the previous section. Finally ao{p,m^) is the cross 
section for the absorption of a dipole at impact parameter p by the nuclear target. This is a non- 
perturbative quantity that normalises the structure function and introduces dependence on the 
mass scale characteristic of the nucleon target. 

So far we have derived the evolution equation that determines the dipole density n(Y, r, p). For 
completeness let us briefly review the calculation of defined in eq. (|J). In momentum space 
and to 0{a s ) the 7* wave function is [fL6|] 



1 roo Ah- 

4 L (z, k) = / — uj+G^k, qh + v e^ L (q) , (47) 

where u and v are the spinors for the outgoing fermions, ej. L {q) is the polarisation vector of 7* and 
is the 7*g<? Green's function 

q) = ieZ « ki-ml + l e ^ {k-qy-ml + l e ' (48) 

with Z q the electric charge of the quark. Performing the numerator algebra and the k~ integration 
we obtain 
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where 



[i 2 = z(l - z)Q 2 + m\ . (51) 



After taking the impact parameter transformation defined by 



r ) = / e ikr 4 0) L (z, k) , (52) 



and squaring as in eq. (Jsj) we find 

2N e 2 Z 2 

4°W) = -^^{[z 2 + {l-z) 2 ]p 2 K l {pr) 2 + m 2 q K {pr) 2 } , 

*i°W) = 8 -^^z 2 (l-z) 2 Q 2 K ( P r) 2 . (53) 

These expressions have also been derived in ref. || from the imaginary part of the 7*5 forward 
amplitude. 



The factorised expression for the structure functions Ft,l, eq. (f46|), allows us to study the 
effects of soft radiation in the small- 2; region beyond the leading power of Q. To this end we shall 
follow the standard lore of renormalon analysis @ . This means that we shall study the singularity 
structure of the Borel image Ft,l{x, Q 2 ; b) defined in the usual way, 

b/a s (Q 2 ) 



F T , L (x,Q z )= dbF(x,Q z ;b)e-»< a ^>. (54) 
Jo 

We emphasise that even in the leading power in Q analysis the structure functions always contain 
a non-perturbative component, denoted by 00 here. However, we are looking for the Borel singu- 
larities that are generated in the subasymptotic Q regime by the components that are calculable in 
perturbation theory, <J?(°) and n in our case. This is how perturbation theory is assumed to signal 
the presence of power corrections. Inspecting eq. ( |46| ) we see that the source of Borel singularities 
is the dipole density n. has no a s dependence and therefore it stays unaffected by the Borel 
transformation. It affects the number density of the produced dipoles N(Y, r, p) by determining 
the distribution of the initial q-q dipole in impact parameter (transverse size) space. <&(°) contains 
infrared regulator set by the scale fj,, eq. (|5l"D. This follows from eqs. d53| ) and the asymptotic 
behaviour of the Bessel--fC functions for fir — > 00, 

Note that the leading term in this asymptotic expansion is independent of the 7* polarisation. 
Since $>(°\z,r) is dominated by values of z away from the end points z = and z = 1, we can 
think of the initial 7* — > qq fluctuation as having size of 0(1/Q). The object that regulates the 
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emission of dipoles of large sizes at the end of the cascade is the dipole-nucleon cross section ao- 
Indeed, dipoles of size p 3> Rn, with Rn = C(l/mjv) the nucleon size, will not couple to the 
target. From this discussion we can formulate the problem of the Borel singularities of the small- x 
structure functions as follows. The Borel image of the structure functions is given by 

Ft,l(x, Q 2 ;b) = f 1 dz f d 2 r ^{{z, r) a dN (Y = ln(z/x), r; b) (56) 

where & d N is a solution to the evolution equation 
8 r°° r b 

w <r dN (Y, r; b) = J dr' J db' £(r, r'; b') d dN (Y, r>; b-b'), (57) 



with kernel given in eq. (21) and boundary condition 

a dN (Y = 0, r) = 5(b) a (r,m N ). (58) 

In the previous section we studied the action of the kernel K. on the test functions (r 2 ) 7 , see 
eqs. (M, H). Then these funct ions where used as a basis for decomposing the dipole density 



in terms of the spectral amplitudes, eq. (^), and we observed that this leads to 7-dependent IR 
renormalons. If this procedure were applicable for constructing the structure functions Ft,l it would 
yield a leading IR renormalon at b/3o = 7, which could be to the left of unity for some < 7 < 1, 
corresponding to power corrections of 0(\/Q p ) with p < 2. This would be in contradiction with 
the standard OPE expansion of the structure functions, which predicts leading power corrections 
of 0(1/Q 2 ), coming from the next-to-leading twist operators. 

In the case at hand such a contradiction does not arise for the following reason. We have noted 
that IR renormalon singularities come from the integration region r' > r in eq. (^), i.e. from 
emission at large impact parameters. This means that in eq. ( |57| ) we are sampling large values of 
r'. Although ao(r' , mjv) cannot be calculated in pQCD it is controlled by a nucleon wavefunction 
which is suppressed at large r', indicating that the probability to find a sufficiently large primary 
dipole inside a nucleon is negligible. We can model o"o (as we shall do in the next section) by 
assuming that it has the same asymptotic impact parameter dependence as that of a virtual photon 
and mass scale p, — > m^. Such a functional form, see eq. (55), cannot be decomposed in the (r 2 ) 7 



basis. We emphasise that this does not mean that IR renormalons do not arise from real emission. 
It means that the leading IR renormalon is anticipated at b@o = 1 leading to 0(1/Q 2 ) corrections 
consistent with the Wilson OPE expectation. 

We conclude this section with the following remarks. The conventional (or Wilson) OPE ex- 
pansion for large Q cannot be used for small- 2; resummation, as we noted in the introduction. 
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Semihard processes involve two large scales and a generalisation of the Wilson OPE is required. 
Such an expansion would not only separate soft from hard particles but also high rapidity from 
low rapidity ones. This type of formalism has been brought to a considerable degree of maturity 
by Balitsky in ref. [17], although it is still far from the point of being a computational algorithm, 



like Wilson's OPE. In the high energy OPE, evolution equations are non-linear beyond the LLA(x) 
approximation, whereas eq. ( p2| ) is linear. This is because our evolution equation does not resum 
the full set of the next-to- leading logarithms of 1/x but only a subset that involves the running of 
the coupling. Hence, our solution is not expected to be unitary but it will parametrise sensitivity 
with respect to low transverse scales. 

5 Numerical results 

In this section we demonstrate by numerical means that the Borel transform of a deeply inelastic 
scattering structure function is regular in the region < bfto < 1 and has a singularity at b/3o = 1. 

In perturbation theory we cannot calculate the cross section o"o of eq. (^) for electron-proton 
scattering. Nevertheless we know that it must be controlled by the nucleon wave function, which 
vanishes rapidly (assumed exponentially) for sufficiently large impact parameter. A reasonable 
model, therefore, is to take the impact factor dependence calculated for the virtual photon, eq. (|53|). 
For the realistic case of D.I.S. this is indeed a model, whereas if we were considering the scattering 
of two virtual photons (onium-onium scattering) then this would be exact, up to a factor arising 
from the integration over the longitudinal momentum fraction, z, in the wavefunction of the target 
photon. We require that the model nucleon wave function be normalisable and this leads us to 
choose the case of a longitudinal virtual photon since for the transverse photon the behaviour of 
the wave function for small impact parameter gives rise to an ultraviolet divergence associated with 
the photon wave function renormalisation. The scale p in eq. (|5^) is set to a typical hadron scale 
of p = m 7v = 1 GeV. To this end we set the Borel transform of cr^(p) in eqs. ( |l4| |58| ) to 

~of N (p-b) = p 2 K (m N p) 2 5(b). (59) 

We consider this modeling of the dipole-nucleon cross section to be appropriate for studying the 
IR effects arising from emission and subsequent interaction of large size dipoles. Note that, as 
stated in the previous section, the large r asymptotics of the 7* wave function, which simulates 
the nucleon target here, is independent of polarisation. So our choice of Kq Bessel function in the 
previous equation is convenient and plausible. 



The form of the input functions prevent an analytic solution to eq. (57) from being found so we 
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perform the integrations over u, r' and the convolution in b' numerically using standard quadrature 
methods. More precisely we expand 

a dN (Y,r;b) = J2^ { dN(r;b), (60) 
so as to bring the evolution equation in the form 

*dN(r;b) = J dr> j* dbX(r, r'; b - b')^ (r';b f ). (61) 

The first convolution of the kernel with the longitudinal input function, c^yv( r '^) exhibits the 
structure we expect from section 4, see fig. (H). We observe that the leftmost singularity in the 
Borel plane appears at b/3o = 1 and we identify this as the leading IR renormalon. There is no 
evidence of any singular behaviour in the region < bfio < 1 which would have been in contradiction 
with the OPE expectation. No singularity emerges as b — ► either, due to the cancellation between 
real and virtual parts of the kernel as discussed in section 3. On the other hand, the expected 
leading singularity at 6/3o = 1 can be seen clearly from the right hand graph of fig. (|l|). For the 
first iteration of the Borel transformed kernel, this singularity occurs as a pole and it is generated 
entirely by the virtual correction contribution to K,. 

There is a distinct difference between the asymptotic behaviour of a^(r;b) and ^^{r;b) as 
fir — > oo, fig. (Q). Whereas <5"d7v( r '^) decays exponentially at large fir, eq. (|55|), the asymptotic 
behaviour of <x^(r; b) is: 

^(r^J^^r) 2 ^- 1 ). (62) 

Examination of the kernel reveals that K, scales as (/xr') 2 ^ 0-1 at large fir' and this is the scaling 
behaviour which dominates in the convolution to determine the behaviour of cr^. For intermediate 
values of fir, the data obtained can be fitted accurately if we include logarithmic correction terms 
and we find 

b) ~ Qtxrf^ ^ («i + a 2 ln( M r) + a 3 ln( M r) 2 ) , (63) 

where oj, a 2 and a% are fit parameters. This fit is exhibited in fig. @ . 

In section 4, we demonstrated that the exponential nature of at large fir was responsible 
for the Borel singularity structure observed in the first convolution. The power like behaviour of 
<5jjl(f*; b) at large fir does not imply that new singularities will appear in the Borel plane to the left 
of 6/?o = 1. It is the action of the convolution in b' in each subsequent iteration which guarantees 
this. For large r' , the convolution integral scales as r'( 2b A'^ 1 ) from the kernel, r' 2 ^ b ~ b )^ 0_1 ) from 
the behaviour of a^(r';b) and the integration introduces a factor of r; the result being that the 
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second convolution also exhibits a power dependence of the form given in eq. (62) for large r. Indeed 
this large r behaviour is sufficient to predict the position of the leading Borel plane singularity for 
subsequent iterations. Since the infrared singularities are determined by the infrared (large r) 
behaviour of the function on which the kernel acts, we can consider 

db' j dr'K{r,r',b-b'){r l2 ) b ^- 1 = -± J db' X {b'Po - 1, b - b') I ^- j (r 2 ) 6 *" 1 . (64) 

Then, from eq. ( p5|) we see that x(6'/?o — 1,6 — b') has a leading singularity at 6/3o = 1, which is 
where the leading singularity will be found for all subsequent iterations. This is demonstrated in 
fig. (|3|). We also note that the convolution in b leads to a linear fall off as 6/?o — ► 0. For further 
iterations of the kernel it is the intermediate 6/?o range which will be of most significance. Although 
the leading singularity occurs again at bfio = 1 as expected from the discussion above, the shape 
of the distribution is broader. This is because for the second and subsequent iterations the leading 
singularity is converted into a cut with branch point at 6/?o = 1 and receives singular contributions 
from both the virtual correction and real emission part of the kernel. 

The generation of each iteration is a numerically intensive task. After generating a sufficient 



number of iterations, we will be in a position to solve eq. (61) and construct the structure functions 
numerically. As explained above, we do not expect the position of the leading singularity to shift. 
However, the behaviour of the wave function in the vicinity of the singularity contains important 
information about the contribution to structure functions from the infrared regions of transverse 
momentum space. We will report on findings for the structure functions in a future publication. 



6 Summary 

Although there exist numerous parametrisations of experimental data for semihard processes that 
include non-perturbative effects in plausible ways, here we have attempted to study in a systematic 
fashion the emergence of non-perturbative corrections as signaled by perturbation theory itself. 
For small- x D.I.S. we have worked within the dipole cascade formalism and introduced the running 
coupling in a self-consistent way. To study the resulting loss of scale invariance we considered the 
dipole evolution equation in Borel (b) space. The variable b provides a measure of deviation from 
the conformal limit 6 = 0, where we have shown that the evolution equation reproduces the well 
known BFKL result. 

We have identified the Borel singularities of the dipole density and have studied how these 
singularities change once the dipole density is convoluted with the dipole-nucleon cross section. In 
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this case the evolution equation generates leading singularity which is a branch cut at bfio = 1. 
This was established by numerical calculation for D.I.S. at small x. Such a singularity indicates the 
presence of 1/Q 2 power corrections to the small- x structure functions coming from the emission of 
dipoles of large transverse size. These correspond to 'ladder gluons' of small transverse momentum 
in the BFKL formalism. Hence our approach is suitable for studying the region of diffusion of 
transverse momentum towards the infrared. Moreover, it would be interesting to see how the 
scale dependence of the interactions in the dipole cascade modify the dependence of the structure 
functions on 1/x. 

Acknowledgement: The authors would like to thank G. Ross for useful discussions. One of us 
(KDA) acknowledges the financial support of PPARC. 
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Figure 1: The Borel structure of a^(r;b), for fir = 1. The graphs show intermediate bflo and 
bPo — > 1 behaviour. 
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Figure 2: The behaviour of a dI ^(r;b), in the intermediate \ir range and at bfio = 0.25, using the 
fit in eq. (pi). 
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